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' We study the ground state properties of trapped atomic condensates with electric field induced 

, dipole-dipole interactions. A rigorous method for constructing the pseudo potential in the spirit 

' of ladder approximation is developed for general non-spherical (polarized) particles interacting 

04 ' anisotropically. We discuss interesting features not previously considered for currently available 

, alkali condensates. In addition to provide a quantitative assessment for controlling atomic interac- 

^ ' tions with electric fields, our investigation may also shed new light into the macroscopic coherence 

I properties of the Bose-Einstein condensation (BEC) of dilute interacting atoms. 
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, The success of atomic Bose-Einstein condensation 
O ' (BEC) has stimulated great interest in the prop- 

O [ erties of trapped quantum gases. In standard treatments 
$_H ■ of interacting quantum gases, realistic inter-atomic po- 
tentials V{R) are replaced by contact forms uoS{R) in 
^ ' the so-called shape independent approximation (SIA) 
^ ; ■ Such an approximation results in tremendous simplifica- 
d . tion. To date, the SIA has worked remarkably well as 
^ ' recent theoretical investigations Q have successfully ac- 
counted for almost all experimental observations 



■ Currently available degenerate quantum gases are cold 
' and dilute, with interactions dominated by low energy bi- 
Q nary collisions. When realistic interatomic potentials are 
I— 1[ assumed to be isotropic and short ranged, i.e. decreasing 
faster than asymptotically for large interatomic 

separations R, the properties of a complete two body 
collision is described by just one atomic parameter: age, 
' the s-wave scattering length. The scattering amplitude 
is isotropic and energy- independent: f{k,k') — — 47rasc 
for collisions involving incident momentum k scattering 
into k' . Effective physical mechanisms exist for control 
of the atom scattering lengths If implemented, 

"Y^i ' these control 'knobs' allow for unprecedented compari- 
' son between theory and experiment over a wide range of 
^ . interaction strength. Indeed, very recently several groups 
I ' have successfully implemented Feshbach resonance ^A\ , 
' O thus enabling a control knob on age through the changing 
of an external magnetic field. Other physical mechanisms 
also exist for modifying atom-atom interactions, e.g. the 
shape resonance due to anisotropic dipole interactions in- 
side an external electric field ]T^ . 

Although fermions with anisotropic interactions are 
' well studied within the context of ^He fluid and in d- 
wave high Tc superconductors, anisotropically interacting 
bosons have not been studied in great detail. In partic- 
ular, we are not aware of any systematic approach for 
constructing an anistropic pseudo potential 1^. 

In this paper, we study the ground state properties of 
trapped condensates with dipole interactions. A rigor- 
ous method is developed for constructing the anisotropic 
pseudo potential that can also be applied to future po- 
lar molecular BEC |14|,p^. This Letter is organized as 
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follows. First we briefly review the SIA pseudo-potential 
approximation. We then construct an analogous effec- 
tive low energy anisotropic pseudo-potential. Numerical 
results are then discussed for *^Rb inside the external 
E-field in the JILA TOP trap. We conclude with a brief 
discussion of prospects for realistic experiments. 

For N trapped spinless bosonic atoms in a potential 
Vt(r), the second quantized Hamiltonian is given by 



2M ^ > ^ 



'l'(r) 



where ^(r) and ^^(r) are atomic (bosonic) annihilation 
and creation fields. The chemical potential fi guarantees 
the atomic number N = J dr^^ {r)'^{f) conservation. 

The bare potential V{R) in (|^) needs to be renor- 
malized for a meaningful perturbation calculation. For 
bosons, the usual treatment is based on field theory and 
is rather involved |p|,p^|-p^. Physically the SIA can be 
viewed as a valid low energy and low density renormaliza- 
tion scheme. The physics involved is rather simple: one 
simply replaces the bare potential V{K) by the pseudo 
potential uo5{R) such that whose first order Born scat- 
tering amplitude reproduces the complete scattering am- 
plitude (— asc). This requires uq — Anh^asc/M. 

When an electric field is introduced along the positive 
z axis, an additional dipole interaction 

Ve{R)= -U2 ^3 , (2) 



appears, where U2 — 4-\/(7r/5) a{0)a* {0)£^ , with a(0) 
being the polarizability, and £ the electric field strength. 
As was shown in Ref. ||T^], this modification results in a 
completely new low-energy scattering amplitude 

/(fc,fc')[^^,^^ = 47r J2 t't'inYLim'rrAk'), (3) 

lm,l' m' 

with t\™ {£) the reduced T-matrix elements. They are 
all energy independent and act as generalized scattering 
lengths. The anisotropic Ve causes the dependence on 
both incident and scattered directions: k and k' = R. 

A general anisotropic pseudo potential can be con- 
structed according to 

V^ffiR) ^ uo5{R) + J2 Ihrnfi^^^. (4) 

l\ >0,mi 

whose first Born amplitude is then given by 

/Bor„(fc,fc') = -{\'Kfa,,YS^(k)Y^S') - ^ E 7hnn(4^)'EE^-" {h.rm)Y:^[k)Yv^,[k'\ (5) 

/imi hn I'rn' 

with Ti;;^\h,m,) = (^)'+''7^f Both 
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Iiy{limi)= {Yi,^,\Yi„n,\Yi,n}, and 

/ 1 

n\= dR-ji{kR)ji,{k'R), 

can be computed analytically |^ . The form in Eq. 

assures all TZ\ tohe k — k' independent [by a change 
of variable to x = kR in the integral] . Putting 

fBorn{k,k') = f{k,k'), (6) 

one can solve for the ji^miiS) as {£) are known nu- 
merically [^,^. This reduces to the linear equations 

- E li^,nA^^)Ti:r' {Km^) EE (7) 

/imi 

for all (/to) and {I'm') with I, V ^ 0, and separately 
asc(^) = ^'^oo(^)- The problem simplifies further for 
Bosons (fermions) as only even (odd) (/, V) terms are 
needed to match. Figure |l| displays result of asc{£) for 
the triplet state of ^^Rb. The Born amplitude for the 
dipole term Ve is 

/B_(fc,fc')=«2-^(47r)2r„2o rl::'Y;^{k)Y,^,{k'), (s) 

Lrn,L m 

with To^o = -0.023508. Tj^' = 7;^"'(2, 0)ITi§ are tab- 
ulated below for small (/, /'). 



TABLE I. Tl^' 
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We found that away from shape resonances, Table 
(|) agrees a few per cent) with the same ratios 
^irn {^)/Jm{^) from the numerical multi-channel calcu- 
lations p2[ . This interesting observation applies for all 
bosonic alkali triplet states we computed: ^Li, '^^^''^K, 
and 85.87j^^^ foj. ^^i^ strength of 3 x 10^ (V/cm) 

|^,|9|. Physically, this implies that effect of Ve is per- 
turbative as £ remains small in atomic units. What is 
remarkable is that Tqq{£) and too('^) ^^^o agree in abso- 
lute values 0]. For ^'^Rb, we found 

U2-^{'i7r)%l° = -1.495 X 10'"t{ao), (9) 

with £ in atomic units (5.142 x 10^ V/cm). qq is the 



Bohr radius. While multi-channel scattering gives |19 

{4Tr)tl"o = -1.512 X 10^"£^{ao). (10) 

The cause of this slight difference (1%) is not entirely 
clear and but is within numerical error. 

We can thus approximate Eq. (^ by keeping only the 
li = 2, mi = term in the sum 

VMR) = uoSiR) - U2Y2oiR)/R^ (11) 

away from the shape resonance. At zero temperature the 
condensate wave function 'il){r,t) — {^{f,t)) then obeys 
the following nonlinear Schrodinger equation 



-|i-V2 + l/,(f)-/i + uo|^(f,i)p 



I ,^Y2o{R) , . ,^|2 

'U2 / dr ^3 mr ,t)| 



^{r,t), (12) 



with ip(f,t) normalized to N. The ground state is found 
by steepest descent through propagation of Eq. (|l^) in 
imaginary time (it). For a cylindrical symmetric trap 
Vt(f) = M{u;'j^x^ +uj]^y'^ +uj'^z'^)/2, the ground state also 
possesses azimuthal symmetry. Therefore the non-local 
term simplifies to 

dr^mp',zT^^- J dz'dpX{.,.;.M{p',z')f, 

with the kernel JC{p, p'; z — z') expressed in terms of the 
standard Elliptical integrals E[.] and K[.]. The kernel is 
divergent at r = r*, so a cut-off radius Rc is chosen such 
that K,{p',p,z' — z) ~ Q whenever |r — r*| < R^. We 
typically R^ ~ 50(ao), much smaller than the grid size, 
to minimize numerical errors. Technical details for nu- 
merical computations and for handling the singular rapid 
variation of the kernel over small length scale will be dis- 
cussed elsewhere Q . 

Figure || presents V(p, z) along p = (a) and z = 
(b) cuts respectively for ^'^Rb (osc = 5.4 nm) at sev- 
eral different £. We note the condensate shrinks radially 
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while stretches along z-axis to minimize the dipole in- 
teraction Ve- The top-right corner inserts shows electric 
field polarized atoms in (radially) repulsive (a) and (lon- 
gitudinally) attractive (b) configurations. An elongated 
condensate along the z-axis reduces the total energy. The 
same mechanism could cause spontaneous alignment of 
polar molecular condensates inside isotropic traps p4[ . 
For better insights we try a variation ansatz 



z) = ^3/4^3/2 exp 



(13) 



with parameters d and k. In dimensionless units for 
length {a± — yJhjMjJZ)-, energy (?iw_L), and A = uj^juj^, 
we obtain 

i^[^TH(l + ^)rf^ + (l + y)^ + ^^^, (14) 

with the effective scattering length a^^ = asc[l — 
&(k)w2/uo], and 



Svr / 2 1 3/*^ tanh ^ \/T~~«? 



h{K) = —^^ -2n' - 1 



3(^2 _ 1) I ^ 



K 



.2 



The &(k) is monotonically decreasing, and bounded be- 
tween b(0) = V5n/3 and 6(oo) = —2^/5tt/3. a^f is shown 
in Fig. ^ as a function of £ for several different values of 
K. For increasing electric field £, variational calculation 
results in decreasing k, eventually k becomes less than 
one, i.e. the condensate changes from oblate (pancake) 
shaped at zero field (for the TOP trap) to prolate (cigar) 
shaped. We also note that 6(k) > for k < 1, therefore 
flgf becomes negative at certain field value £c in the case 
of a positive asc{£ = 0), causing the collapse of the con- 
densate. This is indeed what we found as illustrated in 
Fig. IJ. A detailed discussion of the collapse and other 
interesting features will be given elsewhere pO| . 
We note the energy of dipole alignment 

Ep (27r)l X 10^^ X £^(Hz), (15) 

becomes much larger than the trap depth at the proposed 
£ values for *^Rb. Therefore spatial homogeneity for £{r) 
is required. At £ - 5 x 10^ (V/cm) [£ ~ 10"''] with a 
spatial gradient < 10~'*/cm^, the corresponding force is 
smaller than the magnetic trapping force for typical traps 
at ~ 100 (Hz). For comparison, the magnetic field gradi- 
ent is ^ 10~^/cm'^ inside the Penning trap magnets. Al- 
though the proposed electric field [10^ (V/cm) < £ < 10^ 
(V/cm)] is large, it can be created through careful labo- 
ratory techniques as breakup is fundamentally limited by 
field ionization, which typically occurs at £ > 10'' (V/cm) 
||. Recently a £ field of upto 1.25 x 10^ (V/cmi was 
used successfully to decelerate a molecular beam p3|. 

In conclusion, we have developed a general scheme for 
constructing effective pseudo-potentials for anisotropic 
interactions. Our scheme guarantees that the first 
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order Born scattering amplitude from the pseudo- 
potential reproduces the complete scattering amplitude 
obtained from a multi-channel computation including the 
anisotropic dipole-interaction, thus contains no energy 
dependence at low temperatures of the trapped atomic 
gases ll^. Our scheme is thus more pleasing than the 
standard Skyrme type velocity dependent effective po- 
tentials commonly adopted in nuclear physics ||25| . We 
also presented results for both the electric field modified 
atomic scattering parameters and the induced changes to 
the condensate for ^^Rb in the JILA TOP trap. Our the- 
ory can be directly extended to systems involving mag- 
netic dipole interaction of atoms/molecules in a static 
magnetic trap and systems of trapped molecules with 
permanent electric dipoles For alkali atoms, typ- 

ical magnetic dipole interaction is weak since a Bohr 
magneton [fis — eh/2mc) only corresponds to an elec- 
tric dipole of ~ (l/2a/)(eao) (fine structure constant 
a/ « 1/137), which is equivalent to the induced elec- 
tric dipole at f = 6 X 10^ (V/cm) for ^^Rb. Other 
atoms with larger magnetic dipole moments 26 1 will dis- 
play clearer anisotropic effects. Typical hetero-nuclear 
diatomic molecules have a permanent electric dipole mo- 
ment of ~ (eao), corresponding to an induced moment in 
8^Rb at 8 = 1.6 X lO"^ (V/cm) @. Trapped molecules 
with aligned permanent electric dipoles (by an external 
E-field) would give similar results. However, magnetic 
trapped molecules with unaligned electric dipoles in- 
teracting with the spin axis represents an interesting ex- 
tension that requires further investigation. 

We thank Dr. M. Marinescu for helpful discussions 
during the early stages of this work. This work is sup- 
ported by the U.S. Office of Naval Research grant No. 
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FIG. 1. The field dependent value for Osc- Note the shape 
resonance for £ around 8.3 x lO'' (V/cm). 



FIG. 2. (a) V(P,0) for *''Rb with oj± = {2n)70 (Hz), 
oJz = VS(jJ±, and N = 5000 atoms. Solid, dashed-dot, 
(lashed, and dotted lines are for f = 0, 4.0 X 10^ 5.7 X lO'^ , and 
5.88 X 10'^ (V/cm) respectively, ax is the radial trap width, 
(b) same as in Fig. 2 (a), but for il){Q,z). 



FIG. 3. Typical behavior of the effective scattering length 
a°c • Lines corresponds to k = 5.1, 1.7, 1.02, 0.34, 0.017, in de- 
scending order of al^ . 

FIG. 4. Electric field dependence of the width aspect ratio 
for parameters of Fig. 2. The solid line is the result of our 
variational calculation while circles denot e ex act numerical 
results. The dashed line corresponds to V^\71, the result for 
a non-interacting gas in the TOP trap. 
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